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1. INTRODUCTION 


Because of the high-strength for light-weight capability of carbon 
fiber/epoxy structural materials » a significant increase in their use in 
commercial aircraft structures has been stimulated by the positive effect on 
fuel conservation that the use of such structures would produce. However, 
massive utilization of new materials must be preceded by a thorough examination 
of any potential hazard created by their use. In the case of carbon fiber, 
materials, a potential hazard connected with the release of fibers into the 
atmosphere in an aircraft accident resulting in fire .was identified. 

A necessary part of the hazard evaluation (which has been completed 
and is reported in Ref. 1) is the determination of the detailed characteristics 
of large liquid-fuel pool fires. In a possible post-accident scenario, fibers 
are released through consumption of the epoxy matrix in the lower regions of 
the fire, and are transported up in the fire plume, until they become entrained 
by atmospheric winds. The effective source of fibers from the standpoint of the 
atmospheric dispersion process can thus be significantly above ground level, 
but the altitude to which a given fiber is carried depends on both the fiber 
characteristics and the characteristics of the fire itself. Further, the fibers 
can themselves be consumed within the fire, with the amount consumed dependent 
on the local temperature and fiber residence time in the interior of the fire 
plume. 

To provide the detailed prediction of the characteristics of a large 
liquid-fuel pool fire that is necessary for adequate determination of the 
distribution of fibers that would be entrained by atmospheric winds, a 
mathematical model of a liquid fuel pool fire was developed by SAI. The 
characteristics of the mathematical model were dictated both by the detail 
required to obtain the desired result and by constraints on the time and 
resources available to complete the analysis. Since fiber distribution and 
consumption were necessary results, a relatively simple, one-dimensional fire 
model such as has been reported by Brzustowski (Ref. 2) and by others was not 
suitable. On the other hand, the general solution of the liquid-fuel pool 
fire problem involves the numerical integration of the three-dimensional 
time-dependent, elliptic, Navier-Stokes equations with finite-rate chemical 



kinetics, and the general problem is further complicated by the fact that the 
rate of vaporization of fuel is itself a function of the thermal environment 
including back-radiation from the fire to the fuel pool. 

A tractable mathematical formulation of the pool fire problem can be 
developed if a number of fundamental assumptions are made. In particular, 
if it is assumed that the fuel vaporization rate can be Independently specified, 
(and is constant) and that the fire plume is axisymmetric, the problem can be 
formulated as two-dimensional and steady-state. If it is further assumed that 
over the bulk of the region of interest gradients in the dependent variables 
are much larger in the radial direction than in the axial direction, and the 
static pressure is uniform, the parabolic boundary layer form of the governing 
equations can be used. This assumption is certainly not valid in the region 
near the base of the fire, where the ground plane constraint results in 
entrainment at nearly right angles to the fire plume, but can be expected to 
be reasonably well -satisfied in the buoyancy-dominated upper plume. The 
parabolic plume modeled with the analysis described in this report thus 
represents a simplified analog of a pool fire in a no-wind situation. Although 
simplified and specialized to the no-wind case, the plume model incorporates 
considerable detail, including radial and axial velocity, temperature, and 
species mass fraction profiles, as well as detailed fiber distributions in 
all parts of the fire flowfield. Initial fiber flow (insertion) rates are 
based on empirical data for the consumption of the fiber matrix, and are 
assumed temporally but not necessarily spatially uniform. 

V. 

Within the overall plume model, submodels are required for the turbulent 
mixing process, for gas-phase and fiber consumption chemistry, for particulate 
transport (which includes soot as well as fibers), and for thermal radiation 
from the fire to its surroundings. Each of these submodels is described in 
detail in this report. Briefly, the turbulent mixing process is modeled 
using an eddy viscosity approach; gas-phase chemistry is described with a 
completeTCombustion equilibrium model which assumes different products for 
regions of the plume which are at different local fuel -air ratios; fiber 
(and soot) consumption is treated through the use of a one-step global finite- 
rate chemistry model; and thermal radiation is modeled using a simple optically- 
thin gas approximation. Particulate transport is modeled by discretizing the 



ptrticulate size distribution into a number of classes, each of which is 
treated as a continuun specie. The particles are assumed to be in thermal 
equilibrium with the gas phase, and to move in concert with the mean gas 
phase motion, but can diffuse through the flow at a rate which differs from 
the gas-phase diffusion rate. Buoyancy effects are included in the formulation 
of the gas-phase momentum and thermal energy equations, but effects of buoyancy 
on the turbulent structure of the flow are neglected. These latter effects 
are in any case not well understood. 

Initial results obtained through the use of this model were described in 

an interim technical progress report. Ref. S. This work showed that the fire 

problem is highly nonlinear, so that simultaneous solution of the governing 

differential equations describing the transport of momentum, energy, and gas 

phase and particulate species, including the terms representing interphase 

transport of momentum, energy, and species is necessary to obtain an accurate 

characterization of the flowfield. Fiber consumption within the fire was 

shown to occur, and to be especially significant for small fibers. The initial 

• 

spatial location of the fibers at the point of release was shown to have 
negligible effects on the distribution of fibers in the plume above the flame 
zone, within the limits investigated under the axisymmetric assumption. 

At the time that the work described in Ref. 3 was completed, no large- 
scale pool fire data existed with which to compare and validate the fire model. 
Thus although the results obtained appeared to be qualitatively correct, no 
quantitative agreement of the results with experimental data Could be adduced. 
This, and the large number of assumptions inherent in the model (some of which 
were made necessary by the lack of available experimental data) led to the 
definition, by NASA, of an experimental program to measure the characteristics 
of large scale pool fires. In this experimental program, radial temperature 
and species concentration profiles were measured at a number of axial locations 
within the fire plume, for pool fires of 7.5 and IS meter diameters (25 and 50 
feet), from near ground level to 21 meters (70 feet) above the fire. The 
instrumentation and their locations were laid out based in part on the results 
of the initial fire calculations reported in Ref. 3. 

The results of the NASA experimental program, carried out at the White 
Sands (NM) test facility, showed that the initial pool tire computational results 
were qualitatively correct, but that improvements to the model in several 
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areas were required to obtain better quantitative agreeaent. In particular* 
these results showed that both the turbulent mixing and gas-phase chemistry 
submodels required improvement. Improvement with respect to the overall mixing 
rate in the plume was achieved through modifications to the eddy viscosity 
formulation used: although in general higher-order closures (turbulent kinetic 

energy models) can be expected to provide better predictions of the overall 
mixing rate in a turbulent flow than eddy viscosity fonmlations, the additional 
complexity of such models is not warranted in the present context because of 
the number of other assumptions already present in the analysis. The experi- 
mental results also showed that the chemical kinetic processes which take place 
within the fire plume occur in a considerably different manner than is predic- 
ted by the local complete comlxjstion equilibrium model used in the analysis. 

In part this difference is created by finite-rate chemical kinetics effects* 
but turbulent unmixedness seems also to play a major role. While substantial 
changes to the chemistry submodel in the fire plume model cannot be made without 
introducing significant (and again probably unwarranted) complexity* a simple 
"unmixedness** model was incorporated* substantially improving several aspects 
of the comparison between the model predictions and the experimental data. 

Results of this preliminary work on the improvement of the pool fire model 
were described in Ref. 4. For certain variables* such as temperature* good 
agreement between model predictions and the experimental data was demonstrated. 
However* the overall fuel consumption rate predicted by the model was still 
far greater than that observed in the experiments, and the temperatures predicted 
to exist in the highly fuel-rich portion of the fire were substantially 
higher than those measured. Thus an additional model improvement effort was 
undertaken* focusing on improvement to the prediction of species distributions 
within the fire. It was recognized at the outset of this work that complete 
agreement between the model predictions and the experimental results could not 
in all likelihood be achieved* within the limits of the complete-combustion 
equilibrium chemistry formulation assumed. More detailed chemical kinetic 
models can be postulated* but neither detailed nor global kinetic schemes have 
yet been demonstrated for the very fuel-rich conditions that pertain in the 
lower regions of the plume. Even if such chemistry models existed* their 
incorporation would not be warranted because of the crudity of the treatment 
of the fluid-mechanical processes inherent in the model in this region. 
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For these reasons, the model improvement effort focused on the prediction 
(using an equilibrium formulation) of temperature and species distributions 
in very fuel-rich regions, and the refinement of the unmixedness formulation 
described in Ref. 4. 

All of the work carried out to develop a usable mathematical model of large- 
scale open pool fires in a zero-wind environment is summarized in this report, 
which thus supersedes Refs. 3 and 4. A detailed description of the model is 
given in Section 2, including the basic equation formulation, the particle- 
tracking methodology, and the gas phase and particulate chemistry submodels, 
as well as the turbulent mixing and unmixedness submodels, and the radiation 
submodel. The effects of parameter variation in these submodels are also 
discussed.' This section ccncludes with a discussion of the limitations of 
the parabolic formulation of the pool fire problem: while, as noted earlier, 

much more general formulations of the problem can be defined, these more general 
formulations are quite complex and some of the inherent complexities are also 
discussed. 

The results of this work are described in Section 3 of this report. 

These results include predicted profiles of velocity, temperature, and concentra- 
tions of O 2 and other species within the plume. Computations of fires for 
pool siz'es ranging from 3.75 to 15 meters (12.5 to 50 feet) are presented, 
to provide insight into the expected scaling of pool fire phenomena. Where 
possible, the .predicted profiles are compared with experimental data obtained 
from the NASA/White Sands test program. Section 4 then provides results 
pertaining to predicted carbon fiber distributions and fiber consumption rates 
for the fire sizes computed. 

Overall conclusions of this modeling program are presented in Section 5. 

This work has produced a model of pool fire phenomena which, although approximate, 
is capable of providing considerable detail particularly with respect to expected 
temperature distributions. While the model has been compared with only a 
limited amount of experimental data, these data provide the most complete 
experimental detail on large-scale pool fires that is currently available. 
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NOMENCUTURE 


a 


2 



H 

J 

k 




"j 

N 

P 



Pr 



coefficient in expression for turbulent kinetic energy dissipation 
rate, Eq. 44 

diameter of particulate specie 
turbulent kinetic energy per unit mass 
monochromatic emissive power, Eq. 46 
coefficient in unmixedness model, Eq. 45 
gravitational acceleration 
particulate species mass fraction vector 
'static enthalpy 

vertical distance above fire pool 
total enthalpy 

mechanical equivalent of heat 

constant in Ferri eddy viscosity model, Eq. 10 

constant in Alpinieri viscosity model, Eq. 12 

initial velocity coefficient, Alpinieri viscosity model, Eq. 12 

monochromatic volumetric absorption coefficient, Eq. 46 

equilibrium constant, Boudouard reaction, Eq. 39 

turbulent viscosity length scale 

mass of particulate specie 

particulate number density vector 

pressure 

partial pressure of oxygen 

partial pressure of carbon monoxide 
partial pressure of carbon dioxide 

turbulent Prandtl number 

radiative heat flux 

radiative heat flux at wavelength X 

radial coordi.nate 

particle radius 

Sauter mean radius 

flowfield half-radius, Eq. 11 
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Or QUAUTY 


" 1 

universal gas constant 
particulate radius vector 

Sc 

turbulent Schmidt number 

T 

temperature 

u 

axial mean velocity 

u' 

axial velocity fluctuation 

V 

radial mean velocity 

v» 

radial velocity fluctuation 

W 

molecular weight 


production rate of i^^ gas-phase specie 


reactions 

production of i^^ gas phase specie from 

*F 

W 

Pj 

production of particle class 

X 

axial coordinate 

Subscripts 

c 

centerline 

e 

ambient 

S 

gas -phase 

ij 

species indices 


.th 


particle class 


particulate-phase 


Greek Symbols 

a nolar atom concentration 

8 species mass fraction 

B element mass fraction 

6 ratio of particle to gas phase edd) diffusivity for energy 

scalar turbulent gas phase diffusivity 
turbulent eddy vis.:osity, gas phase 
6 ratio of particle to gas phase species eddy diffusivities 

turbulent dynamic viscosity 
p density 

a ratio of particle f) gas phase eddy viscosity 

turbulent shear stress 
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2. TECHNICAL DISCUSSION 


ORIQINAL PAGE tS 
OP POOR QUALITY 


The modeling of the physical and chemical processes occur ing in a pool 
fire with the presence of carbon fiber particles involves consideration of a 
number of basic physical mechanisms. These mechanisms include diffusion, 
inertia, viscosity, chemistry, chemical kinetics, buoyancy, turbulent transport, 
and the transport of particulate material. Following the modeling philosophy 
outlined in the Introduction, we restrict our attention to the axisymmetric 
(zero-wind) fire problem, with a specified rate of fuel vaporization. With 
this restriction, the computation becomes one of an axisymmetric, buoyant, two- 
phase jet flow. The formulation thus becomes that appropriate for the boundary- 

t 

layer approximation to the equations of motion. 


2.1 DESCRIBING EQUATIONS 

For an axisymmetric, turbulent, two-phase flow including buoyancy effects 
and satisfying the boundary layer apprcximations, the governing equations can 
be written, following Ref. 5, in the following manner: 


Global Continuity 


3rpu 3rpv 
3x 3r 


Axial momentum 

3u 3u 

35e * ^ 


0 


i 1- 

r 3r 




* ] W ^1 


♦ g(pg-p) 


(Ij 

( 2 ) 


in which x and r are the coordinates normal and parallel to the ground plane, 
respectively (see Fig. 1); u and v are the velocity components in the x- and r- 
directions, is the gas phase eddy viscosity (to be described in more detail 
in Section 2.2), p the local density, the mass fraction of all gas-phase 

species in the multicomponent system, the ratio of the particulate phase 

eddy viscosity (diffusion coefficient) for the jth particle-class to the gas 


phase eddy viscosity, 3^ the mass fraction of particulate phase species j, g 
the gravitational constant, and p the ambient density; 
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( 4 ) 


in which 


.th 


W . 3 production of the i gas-phase specie due to homogeneous gas 

g ^ 

phase reactions 

^ i 

production of the i gas-phase specie' from the j particle 
class due to evaporation or heterogeneous reactions. 


Pij 


and 


•F th 

Wp^ = production of the j particle class. 


and in which the 0^ represent the ratio of particle-to-gas phase eddy diffusivi- 
ties of mass and repre ents the gas phase eddy diffusivity of mass; and 


Mixture Energy Equation 


3H 


3H 13 


“ ?3F 


.rpE 


vg 


K * (’« I “j ■ I ■ *) 


/1-E ^ ft 

\ ^ Pk/\ 1 
J "j "pj "pj Pg / ^ **r j 


- i- 2 e. h , e_. 




3B . 
liil _2i 

3 r 


* 6gO - f-) - ^ B j J 

K y J r 





Where 


e. .!Si 
^ 'og 


6. . o. = ^ 

^ ^Tg Si ^ ^Vg 
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in which H represents the total enthalpy of the mixture, h^ the static enthalpy 
of the gas-phase mixture, Pr the Prandtl number, the Schmidt number, the 
static enthalpy of species j, 6^ the ratio of particle-to-gas eddy diffusivity 
for energy, and represents the radiative heat flux from the flow field. 


The mixture energy equation, Eq. (5), may not in general be sufficient 
to define the thermal state of the system. Only in the case of inert particles 
(non-radiating) is there similarity in particle drag and heating processes 
for which "near" dynamic equilibrium would imply "near" thermal equilibrium; 
if "near" thermal equilibrium does not apply, Eq. (5) must be supplemented by 
a particulate phase energy equation. However, in the current calculations, 
the "near" thermal equilibrium assumption has been made. This assumption 
can be justified by the fact that the bulk of the particulate material present 
in the fire plume is soot. Because of the very small size of individual soot 
particles, they rapidly achieve thermal equilibrium with the gas phase even 
though they are responsible for the bulk of the radiative heat loss from the 
fire. The larger classes of fibers, which do not rapidly approach thermal 
equilibrium, constitute only a very small fraction of the flow, in terms of 
both mass and energy loss. It has also been assumed that the values of 6^ and 
appearing in these equations are unity aiwl that Pr = = 1. 

Eqs. (1) - (5) embody several assumptions in the form in which they are 
written. These assumptions include: 


1. The turbulent shear stress can be represented by an dddy viscosity 
model. 

2. The particulate phase species form a dilute, continuum- like suspension 
in the multicomponent gas-phase carrier. 

3. The volume occupied by the particles is negligible. 

4. The particulates are classified in terms of size and composition 
and each class mixes at a rate determined by its ability to respond 
to the dynamics of the gas-phase motion. 

5. Turbulent transport is important only in the direction normal to 
the primary flow direction, i.e. , radially. 


In order to close the system of equations, expressions for the eddy 
viscosity, e^, the species production terms W, and the radiative heat loss i 
must be specified, along with thermal and caloric equation of state. While 
the thermal and caloric equations of state are given by the equations 


11 
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p 

RTZCe./W.) 
i ^ ^ 


and 


H 


g 
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( 6 ) 


(7) 


models for the remaining terms must be supplied. These models must be both 
reasonably accurate and sufficiently simple to allow a tractable computation 
procedure; in view of the other assumptions inherent in the overall pool fire 
model, great sophistication in the formulation of these submodels is not 
warranted. On the other hand, the accuracy of the submodels used must be 
sufficient to obtain a model that will reliably predict the basic details 
of the fire structure over a useful range of fire sizes. Because of their 
importance, each of the submodels will be described in detail in the remainder 
of this section. In addition to the submodels required for the eddy viscosity, 
the species production terms, and the radiative heat loss, the particle tracking 
methodology requires careful development for a situation, such as the carbon 
fiber /pool fire problem, in which the particles are not inert but are themselves 
burning. The techniques developed to account for particulate mixing and 
combustion phenomena will also be described in this section. 


2.2 TURBULENT TRANSPORT MODELING 

As has already been noted, the formulation used in this study includes 
a turbulent eddy viscosity model, i.e., the turbulent shear stress, x,j., is 
written 


T.J. = -pu'v’ 


3u 


( 8 ) 


where is a phenomenological "eddy" viscosity. As originally formulated 
(Ref. 3), the eddy viscosity was defined by an expression postulated by Ferri, 
et al.. Ref. 6, in which 





- ‘’.“J 

e e 


(9) 
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In this expression, tl.e subscript e refers to evaluation in the ambient environ- 
ment and the subscript c to evaluation along the flow field centerline. Since 
for the zero-wind pool fire problem u -0, in this case Eq. (9) reduces to 

“Sg ■ '‘’■ 1/2 ‘’c "c 

where "half -radius,” defined as the distance from the flowfield 

centerline to the radius at which 

pu = 1/2 (p^u^ + p^u^) . 1/2 p^u^ (11) 

Use of Eq. (10) requires specification of the coefficient k: it is character- 

istic of models of the form of Eqs. (9) and (10) (and other models which relate 
the eddy viscosity to mean flow field quantities) that the. coefficient k is 
itself a function of the flow field. However, investigation of the prediction 
of a number of (nonreacting) variable density flowfields using this approach 
showed that reasonably good results can be obtained for a value of k * 0.025 
(Ref. 7) . 

The flow fields utilized in Ref. 7 to investigate the predictive capability 
of the model given by Eq. (10) were neither reacting nor buoyant. Thus it 
is not particularly surprising that while, when compared to the data obtained 
for a IS meter (50 foot) diameter pool fire in the NASA/Nhite Sands experimental 
program, the pool fire calculation incorporating Eq. (10) provided predictions 
in qualitative agreement with the data, quantitative agreement was not achieved. 

This result is shown on Fig. 2, from which it can be seen that the model underpredicts 
the rate of mixing in the early part of the tire development and overpredicts 
the subsequent mixing rate. (The comparison shown in Fig. 2 is with centerline 
temperature, as are most of the comparative results displayed in this section. 
Comparisons of model predictions with other data will be described in Section 3 
of this report.) 

As has already been noted, the coefficient k appearing within the eddy 
viscosity model defined by Eq. (10) is in general a function of the flowfield; 
moreover, the use of a "two stage” viscosity model in which the value of k is 
different in different parts of the flow is a hoary tradition in the application 
of eddy viscosity models to jets. However, physical constraints limit the 
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FIGURE 2. Effect of Variation of Eddy Viscosity Model and Eddy Viscosity 
Coefficients on Centerline Temperature Prediction, 

IS M (50 Foot) Fire. 


range of variation in the constant k: for k 1, the local turbulent shear 

stress level approaches parity with the flow field dynamic head, which is a 
patently implausible situation. Yet, using Eq. (10), values of k approaching 
unity are required to increase the initial mixing rate (given the parabolic 
formulation of the problem) to the levels apparent from the experimental results. 

From visual and photographic observations of the behavior of large pool 
fires it is apparent that these flowfields involve large regions of complex 
recirculating flow near ground level. The visual and photographic observations 
are supported by measurements which show substantial oxygen concentrations 
near the centerline of the fire at heights just above ground level. This 
phenomenon' can clearly not be rigorously modeled under a parabolic (boundary 
layer) approximation. However, an eddy viscosity formulation has been developed 
which provides quantitatively correct mixing rate predictions under just such 
a situation. 

The results of a series of experiments on the mixing of coflowing hydrogen 
(inner) and air (outer) streams have been described by Alpinieri (Ref. 8); 
this work also involved the development of an eddy viscosity formulation which, 
when used with a parabolic mixing model, provided quantitatively correct mixing 
rate predictions for these flow fields. The eddy viscosity model described 
by Alpinieri is substantially different in form from any other locally (mean 
flow) dependent eddy viscosity formulation. An investigation of the details 
of the experiments performed by Alpinieri, described in Ref. 7, showed that 
this flow field exceeds by a substantial margin the criteria established for 
the existence of large-scale recirculation zones within the inner flow. Thus 
the rapid mixing rate observed by Alpinieri (Ref. 8) in the initial region of 
the flowfield is a result not of turbulent diffusion but of large-scale 
convective mixing. The eddy viscosity model described by Alpinieri is then 
not a model for the turbulent diffusivity but (in the initial region of the 
flow) a correlation, which in conjunction with a parabolic analysis, provides 
a substantially better simulation of a flow with large-scale near-field 
recirculation regions. 

The existence of a large scale recirculation region in the nearfield is 
a common feature of the experiments described by Alpinieri (Ref. 8) and the 
observed behavior of a pool fire. Thus the original eddy viscosity model 
described by Eq. (11) was replaced by the expression 
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where k has been used to replace the tern u in Alpinierl's expression, and 

V 0 

taken to be 1.0. k is an adjustable eddy viscosity coefficient; as was the 

case for the viscosity model given by Eq. (11), the values of this coefficient 

nay be different in different regions of the flow. For the particular case 

of the pool fire, for which the initial velocity u^ is of the order of 0.025 

m/sec (0.08 ft/sec), the bracketed velocity term in Eq. (12) is initially of 

the order of 1/u^, compared to a direct u^ proportionality expressed by Eq. (11). 

As the computation continues, u , driven by buoyancy effects, increases to 

c 

^ 20-30 m/sec (66-98 ft/sec), so that the bracketed term becomes of the order of 
1.5 u^. Thus the effect of the use of this model is to provide a substantially 
greater eddy viscosity level early in the fire computation than would be 
achieved by the more conventional approach embodied in Eq. (11). 

Results of the use of the Alpinieri formulation, Eq. (12) are compared 
with the results obtained using Eq. (11) and the experimental data for the 
15 m (50 ft) fire on Fig. 2. Computations with the Alpinieri model using two 
sets of viscosity coefficients are shown. It can be seen from this comparison 
that both the early (low altitude, h) and late (high altitude) results are 
considerably improved through the use of the Alpinieri model; the changeover 
between coefficient values in each of the Alpinieri model computations occurs 
at h « r^ (7.5 m). This improvement in the predictions for the pool fire 
is not an endorsement of the use of the Alpinieri approach in'general, but 
these results indicate that strictly from the standpoint of the use of the 
eddy viscosity model in conjunction with a parabolic analysis as a correlative 
tool, an improved correlation can be obtained. 

The results shown in Fig. 2 also show that the overall centerline tempera- 
ture profile prediction is poor in the region of maximum temperature. Tempera- 
tures near to the adiabatic flame temperature are predicted to occur along 
the flow field centerline, but nowhere in the flame did measured temperatures 
approach this limit. Since the temperatures were measured using unshielded 
thermocouples, some error is inherent in the experimental results, but this 
error could not explain the 700-800"K discrepancy observed. Further, gas 
sample data indicated considerable unburned hydrocarbons in the region where 
maximum gas phase temperatures are predicted. Both of these observations 
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indicate that the effects of chemical kinetics and a large-scale unaixedness 
are of importance in the pool flame development process. 

2.3 CHEMISTRY MODELING 

Because the objective of the pool fire modeling work described in this 
report was to obtain a computation of the structural details of the fire that 
are important to the prediction of carbon fiber transport and consumption, the 
basic model formulation requires an accurate definition of the thermal field 
rather than the details of the chemical species. Thus the selection of a 
chemical system in terms of the number and kind of species considered can be 
somewhat arbitrary. 

In a hydrocarbon- air system the principal products of combustion are H 2 O 
and CO^; for fuel-lean conditions O 2 will remain while in the fuel-rich regime 
unburned fuel will be found with the products of combustion. Dissociation of 
the major species at high temperatures results in the presence of H, 0, OH, 
and CO in the fuel lean region, while under fuel rich conditions even more 
species, including a variety of hydrocarbon fragments and soot, will appear. 

At temperatures above about 2500*K the dissociation of nitrogen will begin 
to become significant. 

All of these species can be included in an equilibrium chemistry computa- 
tion and many of them have also been included in finite rate chemistry models 
(although considerable work is still required to establish usable finite 
rate chemistry models for complex fuels under fuel-rich conditions). However, 
both full -equilibrium and finite-rate chemical kinetics models require 
considerable computational time and involve a level of computational complexity 
that is not warranted given the objectives of the pool fire calculation. 

Instead, as cited above, the principal quantity of interest in the pool fire 
calculation is the temperature, and it has been demonstrated by Edelman et al. 
(Ref. 9) that accurate temperature predictions, for temperatures s 2S00*K, can be 
obtained using a simple quasi-complete equilibrium formulation. This is 
demonstrated by the results shown on Fig. 3, from Ref. [9], in which the 
temperatures predicted by the quasi-complete equilibrium model are compared 
with those obtained with a full equilibrium formulation, for a propane-air 
system. Good agreement can be seen except near stoichiometric and at the 
higher equivalence ratios. 
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FIGURE 3. Comparison of Quasi^Complete Combustion Model with 
Equilibrium Chemistry for a Propane*Air System. 


In the form used in the pool fire model, the spectrum of equivalence 
ratios encountered in the fire is divided into distinct regimes, defined by 
the fuel/oxygen ratio; for fuel rich regions, a temperature limit is also 
used to define regions in which soot is produced. For the fuel -lean regime, 
complete oxidation of the fuel to the products CO 2 &nd H^O is assumed to 
occur. The upper limit for this regime is the stoichiometric point where the 
atom balance is given by 


■ *0 

in which the S's are the element molar concentrations. In terms of element 
mass fractions, this relationship can be written 





( 14 ) 
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Accordingly, the lean regine is defined within the linitt 



> 0 


and the specie mass fractions are given by 
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In the fuel rich region, pure fuel remains if the local temperature is below 
7S0*K; for temperatures above 750*K, soot can be formed. The fuel -rich region 
is itself divided into two zones. For fuel-oxygen ratios just above stoichio- 
metric, the fuel is assumed to be oxidized to CO and COj* while the molar 
concentration of water remains essentially constant. Thus, for 


"c "o-*c 

the species mass fractions are given by 
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a further subdivision is defined depending on the temperature and the carbon/ 
oxygen atom mole ratio. Thus, for T < 750* K and ^ , the products 

are assumed to be COj, and fuel so that ^*0 



(24) 

(25) 
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while for T < 750*K and products are assumed to be Hj, CO, 

and C H , and 
n m 
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If the temperature is above 750*1C, and 8 q/Wq > it is assumed that the 

fuel decomposes and partially oxidizes to CO, H 2 , and H 2 O, so that 
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vhile, if T > 750*K and > Iq/Wq, soot is formed: 

Wp . 

®C(s) “ ■ Wq ^0 

'co ■ 


Since nitrogen is assumed not to react, it satisfies the relationship 


(S2) 


(33) 

(34) 


(35) 


•nj • 

throughout the flow. 

Note that when solid carbon (soot) is formed, it becomes a particulate 
phase specie and thus enters the particle tracking methodology to be described 
subsequently in this section. It may also be consumed in its subsequent 
passage through the fire, producing CO 2 ; the rate fo** this conversion is 
given by the expression obtained by Lee, Thring and Beer [10]' 


dm 


«d 


2 dt 


4 ”2 

1.085 X 10 — exp(-39,300/RT) 
T 


(37) 


in which d^ and m^ are the mean diameter and mass of the j 


th 


particle class. 


and T are the local oxygen partial pressure and temperature, in atmospheres 
ani degrees Kelvin, respectively, and R is the universal gas constant in units 
of cal/g-mole *K. 


Given this chemistry formulation, the computation proceeds in the 
following manner: Transport equations are written for the four gas-phase 

elements H, 0, C, and N, and for each of ten particulate-phase "species," 
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which constitute two classes of soot particles and eight of carbon fibers. 
Solution of these equations, along with the momentum and energy equations 
yields the gas-phase elements at each point, and the chemistry model given by 
Eqs. (19) -(36) yields the species taken to exist at a given point in the 
flow. With the gas-phase species defined along with the local enthalpy level, 
a temperature is obtained using Eq. (7). If soot is formed, the appropriate 
mass fraction is deleted from the gas-phase carbon element and added to the 
appropriate particle class. Under a thermal equilibrium assumption the parti- 
culate species are taken to be at the same temperature as the gas phase, so 
that this temperature is used to compute the particulate consumption rate for 
each species and the thermal radiation from the particulate species, both of 
which will be discussed further in subsequent sections of this report. The 
results of the consumption process represent a source terra for the gas-phase 
carbon element for the next step in the computation. 

Although the chemistry modeling approach outlined in this section was 
reasonably successful in the flame computations described in Ref. 9, the 
results displayed in Fig 2 indicated several unsatisfactory features of the 
model. One of these is the underpr edict ion of the temperatures observed at 
the centerline in the early part of the fire, and the second is the strong 
overprediction of peak centerline temperature. The underprediction of the 
temperature in the early part of the fire can be traced in part at least to the 
chemistry model itself: as seen on Fig. 3, for equivalence ratios greater than 

V. 

about three the basic model of Ref, 9 strongly underpredicts the temperature 
for a propane-air mixture. Both measurements and computational results show 
very fuel-rich equivalence ratios at the centerline for h < 4m, with measured 
equivalence ratios of 10-20 being observed. On the other hand, the overpredicted 
peak temperature is due to the prediction of stoichiometric conditions along 
the fire centerline; with the basic chemistry model temperatures close to 
the adiabatic flame temperature will be predicted at ^ a: i. Since temperatures 
this high were never measured within the pool fire, it is evident that other 
factors were involved in the experiments which resulted in the temperature 
levels observed. 

Thus the comparison of computational results with experimental data such 
as depicted on Fig. 2 dictated further modifications to the basic chemistry 
model. These modifications took two forms. In view of the results shown on 
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Fig. 3, and of other computations, the basic chemistry formulation was modified 
to incorporate a partial -equilibrium expression for the species CO, C02* 

C in very fuel-rich regions. Examination of these and other results and observa- 
tion of the behavior of the pool fire plume during the experimental program 
indicated that there were large-scale unmixedness effects, i.e. , that 
although on the average local fuel -air ratios near stoichiometric were seen, 
this average was made up of the sequential passage of fuel-rich and fuel-lean 
elements, so that the average fuel-air ratio did not adequately describe the 
actual flow. To incorporate this phenomenon, an unmixedness model was 
developed and applied to the calculation. 

2.3.1 Partial Equilibrium for Fuel-Rich Regions 

The chemistry model modification introduced in the fuel-rich region, for 
Ec/Bq ^ ^ ^ 750*K involves the use of the so-called Boudonard reaction 

(Ref. 11) 


CO t y CO 2 ♦ ^ C(S) 


(38) 


as a partial equilibrium relationship. The equilibrium constant for the 
Boudouard reaction is 




and values 


of K‘ 


are quoted in Ref. 


11 and given below: 


T*K 


T*K 


T*K 


300 

1.63x10"^^ 

900 

0.195 

1750 

1.06xl0‘ 

500 

1.82x10"® 

1000 

1.91 

2000 

4.21x10' 

700 

2.71x10"^ 

1400 

628 

3000 

2.8x10^ 


(39) 


For application in the fire code, these equilibrium constant values were 
curve fit by the expression 


In K® » 20.649 - 2. 0002x1 0^/T 

P 


(40) 
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and the soot concentration was conputed using this equilibrium constant rather 
than through an arbitrary specification of the levels of C0« CO 2 and soot. 

Note that Eq. (40) requires the local temperature: since this is a function, 

given the enthalpy, of the species distribution, the Boudard equilibrium is 
computed iteratively in the fire model. 


2.3.2 Turbulent Unmixedness Modeling 

In the simple approach adopted for this effort, it is assumed that because 
of unmixedness a portion of the total oxygen and fuel at each point in the 
flow is unavailable for reaction, and that the fraction unavailable can be 
realated to the local velocity fluctuation level. The proportion of fuel 
unavailable is taken to be related to the unavailable oxygen by the local 
stoichiometry. Thus, although local 4>=1 will still be observed, the temperature 
computed for this point will be modified by the dilution of the products of 
combustion by the amount of oxygen and fuel deemed to remain unburned. 

Therefore, it is assumed that 




t 


where "u" and *'t'* refer to "unburned" and "total" respectively, and 


(41) 
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Then, assuming that 


Hi 

u 


is related to the local turbulent kinetic energy level. 


Hi 

u 



(43) 


and noting that, consistent with the use of a locally dependent eddy viscosity 
model, the turbulent energy production can be equated to the turbulence energy 
dissipation: 


P 




( 44 ) 


a relationship for fr can be obtained: 
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where is the local turbulent viscosity, a length scale taken to be that 
used to obtain the local turbulent viscosity, p is the local density, and f a 
coefficient which is to be evaluated through comparison of predictions with 
data. 

2.3^3 Effect of Chahistry Model Modifications 

Of the two model modifications discussed in the preceding sections, the 
unmixedness formulation has by far the greatest effect. This is because the 
Boudouard equilibrium model is effective in only a very- limited region of the 
fire plume', near ground level. Because of the massive entrainment of air near 
the ground level in the calculated plume, the predicted fuel -air equivalence 
ratio along the axis quickly drops to a level at which the Boudonard equilibrium 
relationship provides only a marginal difference in calculated temperature 
from that obtained using the basic chemistry model. 

On the other hand, the unmixedness formulation has a strong impact on 
the predicted plume centerline temperature development, as is shown in Fig. 4. 
Results for two levels of unmixedness are shown on this figure, in comparison 
with the previously presented results for no unmixedness and with the experi- 
mental data for the 15 meter (50 foot) diameter fire. The major impact of the 
unmixednesness model is to drastically reduce the predicted peak temperature 
along the centerline: in the calculations shown on Fig. 4 this reduction 

amounts to more than 700*K for the model with f = 0.75 (see Eq. 45) compared 
to without the unmixedness correction (i.e., f=0) . In all computations shown 
on Fig. 4, the same eddy viscosity model and value of the emissivity in the 
radiation model (Section 2.4) has been used. Noted that the incorporation 

i 

of the unmixedness effect also changes the shape of the predicted temperature 
curve, particularly at heights below about four meters. This is because the low 
temperatures generated though the incorporation of unmixedness alter the 
bouyancy and thus alter the overall plume flowfield development. The initial 
entrainment of '’mbient air is driven by the assumed initial velocity of the 
fuel, and is constant for all of the cases shown, but from this point on the 
mixing process is driven by the buoyancy term which depends directly ^n the 
predicted temperature level. Since this in turn affects the local velocity 
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FIGURE 4. Effect of Incorporation of Unmixedness Model On 
Predicted Centerline Temperature Distribution 
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gradient, and thus, through Eq. (45) the computed unmixedness, the effect 
of unmixedness on the plume centerline temperature is not a linear function 
of the degree of unmixedness, as these results clearly show. 


2.4 PLUME RADIATION MODEL 

Radiation from the fire plume to the surroundings can be a significant 
source of energy loss from the plume. In order to account for this loss, 
represented by the terra in Eq. (5), the radiation from the plume has been 
computed for the optically thin limit. 

In the optically thin limit, the monochromatic radiant energy loss term 

can be written 

; 


dr 




( 46 ) 


(Ref. 12) in which e^^^ is given by Planck’s law and is the monochromatic 
volumetric absorption coefficient. The term e^^^, which is a function of the 
local temperature, has been written as a function of radius since in the plume 
T = T(r). 

For the present work, Eq. (46) has been used in the form of an average 
(over wave number) of the expression given with k = 0.25. Since the parti- 
culate phase carbon is the most important radiation contributor in the fire, 
gas phase emission has been ignored. Approaches similar to that outlined here 
have been successfully applied to the prediction of carbon particle radiation 
from rocket plumes (Ref. 13). 

Although this approach is simple to apply, the optically thin limit is 
not a good approximation to the behavior of the radiation flux in a dense, 
sooty plume: under these circumstances the radiative transport involves 

the absorption and re-emission of radiation from soot particles all along the 
lines of sight within the plume. However, the optically thin model can be 
used to gain some insights into the effects of radiation by variation in the 
values assumed for the soot particle emissivity, where lower values of emis- 
sivity account to some extent for the absorbing/emitting character of the actual 
fire plume. The results of several computations, incorporating the unmixedness 
correction described in Section 2.3, with f = 0.5, and using the Alpinieri 


eddy viscosity model, Section 2.2, but varying the assumed particle-phase 
emissivity from 0 to 0.9 are shown on Pig. S. These results indicate that 
neglecting radiation effects is not a feasible approach to the computation of 
the plume development, and that for all other effects constant, the predicted 
plume development is a strong function of the radiant energy loss. It is also 
of interest to note that at an emissivity of O.S, the total radiant energy loss 
is about 20% of the total energy produced by the consumption of the fuel, which 
is essentially the proportion of radiant energy/total energy observed in 
laboratory scale pool fires. 

2.5 PARTICULATE TRANSPORT MODELING 

i 

Since the basic objective of the development of the fire plume model 
described in this report was to provide a means for the computation of the 
transport and consumption of graphite fibers released during a pool fire, 
a particulate tracking methodology is a required feature. This particle 
tracking model is also used to compute the transport of soot particles within 
the fire plume. 

The combined effects of continuous reaction of both fibers and soot and 
the non-uniform velocity and temperature fields within the fire plume provide 
conditions resulting in a spectrum of particle sizes at every point in the flow. 
To completely specify the particle reaction rate and the production of gas 
phase species from the consumption of the soot and fibers, the local particle 
size and number density distribution must be known. This in turn requires the 
development of a methodology for tracking the particle history throughout 
the flow. 

The track the particle history through thd flow, the total amount of 
particulate material present at every point in the field is divided into a set 
of categories. These categories can be based on any parameter which can be 
used to differentiate particles; for the fire application particle size 
(diameter for the assumed spherical soot particles and length for the fibers) 
is the relevant parameter. Thus categories are set up by dividing the total 
range of particle sizes into classes defined by a minimum and a maximum particle 
size. 


28 



NE TEMPE 






ORfQfNAL PAGE IS 
OF POOR QUALirr 


Consider, for sinplicity, spherical particles. Within each site range, 
all particles are combined to form an average with a volume to surface area 
ratio equal to that of the summation over all contributions. This equivalent 
particle size is the Sauter mean value and is expressed as 


r ' ■ 
s 


Enirf^/ 


( 47 ) 


where n^ particles have radius r^^. The mass fraction of each category is . 
equal to the sum of the mass fractions of each contribution. With the average 
size of category particles and the total category mass fraction known, the 
appropriate number of equivalent particles is established for every class. 

This representation provides an exact equivalence for the subsequent rate 
process provided the rate of change of mass depends directly upon the particle 
surface area. 

The particle distribution array in terms of size, mass fraction, and 
number in each category at every grip point can only be altered by changes in 
size of the particle due to chemical reaction and by diffusion of particles from 
the neighboring flow. The radius of the particle changes in accordance with 
the particle consumption expression as the flowfield develops. Thus, after 
every step, the new particle size is giver, by 


r ■ r^ ♦ 
o 



( 48 ) 


After every step the updated radius array must be reclassified in order to 
account for the transition of particles from one category to another. Since 
it is possible for more than one contribution to be made to a particular 
category, the total category mass must then be reaveraged using the Sauter 
criteria, to establish the proper radius and number for the group. The particle 
array in terms of mass fraction, size and number density resulting from the 
kinetics of the combustion process will be defined as: 


G(J,I),^ 

R(J,I)^ 

N(J.I)^ 


at grid point 
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Each category in the particle array is altered by the diffusion of 
particles of the same category from neighboring upstream locations. Subdividing 
the particle mass fraction into each separate category and treating the 
category mass fractions as separate species, the diffusion equation can be 
written in finite difference form in von Mises coordinates as 



where the n,m subscripts refer to the grip point locations shown in Fig. 6. 

The second term in Eq. (49) represents the change in due to the 
kinetics process as the flow moves downstream from n to n*«-l. The sum of the 
first two terms is the total mass fraction of category particles present 
downstream due to kinetic processes, and the remaining term represents the 
contribution to the mass fraction from diffusion. Noting Fig. 6, the diffusional 
term can be split into the contributions from each of the three upstream points: 
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(50) 

(51) 

(52) 


For each of the categories present at a given grid point, we may express the 
category mass fraction as 

G(J,I)n+l,m » |ax + AG^ ♦ AG2 ♦ AGj (53) 

where the kinetics contribution is evaluated as G(J,I)^ with its associated 
R(J,I)j^ and N(J,I)j^. Knowing the category sires associated with each upstream 
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FIGURE 6. Finite*Difference Grid System. 


contribution the Sauter criteria (or other appropriate characterization cri- 
teria) can be used to define the new category radius: 


R(J,I)n+l,m » 


G(*I»I)k ♦ AGj ♦ AGj dGj 

... dG, dG. dG- 

^ (d I ^ ) h 11 3 

R(J,I)k R(J,I)n,m*l R(J,l)n,m R(J,I)n,m-l 


(S4) 


Knowing the new mass fraction and size, the number density of particles is 
determined from: 


N(J,I)n+l,m 


pn»l ,mG(J,I)n-fl,m 
4/37rpj^R(J,I)n+l ,m 


(55) 


2.5.1 Effect of Source Location 


A variable of interest with respect to the distribution and consumption of 
fibers in the fire is the effect of the ground level location of the fibers 
within the pool. To investigate the effects of source location, the nonuniform. 



but axisymmctric fiber source locations shown in Fig. 7 were used as initial 
conduions for a 30 meter fire calculation. In all cases the total fiber 
flow rate (10 kg/min) was the same, but for one computation the fiber source 
was assumed to be a centrally-located disc, while in the other case an annular 
ring source was used.* 

The results of these computations are shown in Figs. 8(a) - 8(c). compared 
to the results of a uniform fiber distribution calculation. Although some 
differences in detail are apparent in these figures, it is clear that overall 
the differences are small and would be insignificant in the computation of 
the subsequent atmospheric dispersion process. 

i 

Some insight into the reasons for the insensitivity of the results to 
the initial particle source location can be gained from the isotherm and 
stream function profiles shown in Fig. 9. The isotherm profiles shown a 
strong "necking down" in the early part of the calculation, the reasons for 
which are made apparent in the stream function plot. Since the stream function 
is proportional to the square root of the mass flow in the plume, it can be 
seen that the vast bulk of the flow at the top of the plume results from the 
massive entrainment of air generated by the large buoyancy-driven accelerations 
in the fire. Indeed, the « 1 stream function, which represents the mass 
that originates within the pool, is very quickly driven to the centerline of 
the flowfield by the air entrainment. Thus material originating in the pool, 
whatever its location, is quickly driven in towards the centerline, limiting 
the overall effects of source location. The differences observed between 
source location results in Figs. 8(a) - 8(c) are primarily a result of 
different turbulent diffusion rates developed as a result of the different 
initial distributions. 


* These computations, originally described in Ref. 3, utilize a fire plume 
model which does not incorporate the unmixedness effects described in Section 
2.2. However, differences in fire plume structure observed with the incorpora- 
tion of unmixedness would not materially change the results described in this 
subsection. 
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FIGURE 7. Fiber Source Location* 
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FIGURE 8(a). Effect of Source Location on Distribution of Particulates 
at End of Fire Calculation - Single Fibers 



FIGURE 8(b). Effect of Source Location on Distribution of Particulates 
at End of Fire Calculation • Single Fibers 
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FIGURE 8(c), Effect of Source Locatior. on Distribution of Particulates 
at End of Fire Calculation - Fiber Clusters 


2.6 FIBER CONSUMPTION MODEL 


The Lee, Thring and Beer model (Ref. 10) for the consumption of soot 
particles was utilized in the fire modeling described in Ref. 3 to predict 
fiber consumption rates. In this work, the fiber consumption model has been 
modified to reflect the results of work done at NASA-Langley on individual fiber 
consumption rates as a function of temperature and oxygen partial pressure. 
However, only one set of data, for a temperature range of 450*C to 700*C 
and an oxygen partial pressure of 0.0816 atm was suitable for the development 
of a fiber oxidation rate expression. Thus, while a modified rate expression 
was developed from these data as shown on Fig. 10, it can only be considered 
to be preliminary as the possibly substantial effect of oxygen partial pressure 
on the rate has not been determined. Further, the experiments involve a 
variable particle heat up time during which some consumption is occurring. 

This variable heat up time cannot be incorporated into the development of 
the model for particle consumption. Finally, possible physical differences 
(e.g., surface coatings) between the test particles and those released in a fire 
have not been investigated. All of these uncertainities can lead to substan- 
tial differences between predicted and measured particle consumption rates in 
a fire environment. 
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FIGURE 9. Isotherm and Stream Function Profiles for 30m Fire 
2.7 FIBER THERMAL NONEQUILIBRIUM 

The computations reported in Ref. 3 were all carried out under the 
assumption of fiber thermal equilibrium, that is, that the fibers are 
everywhere the same temperature as the gas phase. This assumption can be 
relaxed, so that the heating and subsequent cooling of the fibers in response 
to the gas phase temperature change is computed. This is done by writing a 
particle-phase energy transport equation for each particle class, which 
includes as source terms expressions representing the heating of the particle 
by conduction from the gas phase, energy loss from the particle caused by 
radiation, and heating of the particle caused by surface combustion, all of 
the energy of which is assumed to be absorbed by the particle. 


EXP. DATA® Pq 2 = 0.0816 atm; Dl AM » 7 in) 
5.83 X 10^ P„2 

OXID. RATE = EXP. (-23,500/RT) 

T^ 

EXP. DATA 

PREDICTION 





hours 


FIGURE 10. Comparison of Modified Lee, Thring, and 
Beer Model With Fiber Consumption Data. 


To assess the need for the incorporation of a thermal nonequilibrium 
model, a number of one-dimensional single fiber computations were carried out 
which are of interest with regard to the response of the fibers to temperature 
changes and the effects of surface combustion on particle temperature. Fig. 11 
depicts the response of a single fiber to gas-phase temperature change consid- 
ering only conduction heat transfer (which is, however, the largest single 
contribution by a substantial amount). For this calculation, the initial gas- 
phase temperature (T^^) was 1500*K, and the initial fiber temperature (Tp^) was 
300*K. Under conduction-only conditions, and assuming that the particle spe- 
cific heat and gas conductivity are constants, an analytical solution can be 
obtained; in the numerical solution both conductivity, X, and heat capacity, 
Cpp, are functions of the gas phase and particle temperatures, respectively. 
Fig. 11 shows that the 1/e response time of the fibe^^ is about 0.10 seconds. 
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Within the fire, there may be regions of locally high (i.e., near stoichiomet- 
ric) temperatures on scales of the order of 0.01 dp, where dp is the local fire 
diameter. If such high temperature regions exist, and a local fire diameter of 
the order of one-half the pool diameter is assumed, and it is further assumed 
that these eddies have a convection velocity relative to the mean velocity of 
the order of O.IU, then a "hot" eddy lifetime is of the order of 0.03 seconds. 
Comparison of this lifetime with the single-fiber response time of 0.10 seconds 
indicates that the fibers are essentially unaffected by local temperature fluc- 
tuations. 



FIGURE 11. Response of a Single Fiber to Gas-Phase Temperature, 

Tp^-Tg^ = 1200*K, Conduction Only. 

Additional one-dimensional calculation results including the effects of 
thermal radiation and particle consumption are shown in Fig. 12. As can be 
seen from this figure, the inclusion of radiation has a negligible effect on 
the temperature history; and the incorporation of pa-^ticle combustion also 
leads, for these conditions, to negligible chanr^s relative to the conduction 
only solution. 


TPo-Tgo»1200“K 







sec 

FIGURE 12. Effects of Radiation and Particle Combustion 
on Fiber Temperature History. 

2.8 COMMENTS ON PARABOLIC MODELING OF LARGE POOL FIRES 

It was pointed out in the introduction to this report that the modeling 
of the open pool fire problem using the boundary layer form of the equations 
of motion. This approach was nevertheless adopted in order to use existing, 
well-developed solution techniques, and because the information desired from 
the model comprised fiber flux and distribution results in the region of the 
fire in which the parabolic approximation is reasonable. The applicability 
of the parabolic formulation and be assessed at least in part by reference 
the stream function plots shown on Fig. 9. Under the parabolic assumption, 
gradients in the dependent variables in the axial direction (i.e., vertically 
along the plume axis) are assumed negligible compared to gradients in the 
radial direction horizontally). The stream function contours shown on 

Fig. 9 demonstrate that near ground level, this assumption is not met, and 
axial gradients can be expected to be the same order of magnitude as (if not 
larger than) the radial gradients. 
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The parabolic approach is thus an approximation to a complete fire model, 
and the errors inherent in the approximation cannot be evaluated without a 
more exact model for comparison. It is, however, clear that the major inaccura- 
cies involved in a parabolic formulation of the pool fire problem are found 
at the lower levels of the fire, and it is in just this region that the major 
discrepancies between the computational results and the experimental data 
are observed. An "exact" formulation of the pool fire problem would involve 
the solution of the full Navier-Stokes equations, including axial diffusion 
and radial pressure gradients, but there are several obstacles that would have 
to be overcome to develop the more exact model . 

Elliptic computations, such as an "exact" formulation of the pool fire 
problem, require the definition of boundary conditions on all of the boundaries 
of the problem. For the zero-wind pool fire, this means that boundary condi- 
tions must be defined a sufficient radial distance from the fire location so 
that the influx (entrainment) specification does not affect the overall 
solution: mathematically the proper boundary conditions are defined at 

infinity. Further, the actual fuel vaporization rate is controlled by the heat 
radiated from the fire to the liquid pool, so that the evaporation rate 
boundary condition is both time-dependent and ill -posed in the sense that it 
depends on the solution for the overall fire structure (which itself depends 
on the boundary conditions applied). The entrainment boundary conditions are 
more easily specified in the case of a fire in a crosswind, for in this case, 
the mean velocity profile approaching the fire can be defined directly, so 
that although the crosswind fire problem is three-dimensional whereas the zero 
wind fire is axisymraetric, the former problem is more easily posed. However, 
the problem of specifying the boundary condition at the fuel surface still 
remains: defining a constant vaporization rate, as was done in the parabolic 
computations reported herein, is at best a crude approximation. 

Thus, although a parabolic pool fire model is a rough approximation 
to reality for the zero -wind case, the development of a more exact approach 
is itself problematical and involves a several orders-of-magnitude greater 
effort, both in development and in solution time. The approximations involved 
in the parabolic formulation improve progressively as height above ground 
increases, so that if what is desired is fire structure and transport data 
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three or more pool radii above the fire surface, the parabolic model should 
suffice. This will be demonstrated, at least in part, by the comparisons of 
model predictions with experimental data to be described in the next section 
of this report. 
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5. COMPARISON OF MODEL RESULTS WITH POOL FIRE DATA 


On the basis of the centerline temperature comparisons shown in Section 2 
of this report, the Alpinieri eddy viscosity approach, with coefficients of 
0.07/0.005, changing at one radius above the liquid pool, using the unmixedness 
model with the factor f « 0.75, and the thin gas radiation approximation with 
an emissivity of 0.50, was selected for detailed comparison with the experimental 
data obtained from the NASA-White Sands test program. It is clear from the 
comparisons shown in Section 2 that, because of the close coupling of each of 
these parameters in the overall pool fire model, large variations in predicted 
results can be obtained with variations in these parameters. Thus the set of 
parameters selected for the comparisons to follow do not represent an "optimum,'* 
but rather represent physically reasonable values which provide a reasonably 
good prediction of the centerline temperature profile observed during the 
experimental program. 

\ 

3.1 RADIAL TEMPERATURE PROFILE PREDICTIONS 

Model predictions of the radial temperature profiles for each of the 
vertical stations at which temperature data were obtained are shown in compari- 
son with the data in Figs. 13-17. At low altitudes in the fire, peak tempera- 
tures are strongly overpredicted by the model, as evidenced by the 1.43 meter 
results shown in Fig. 13, and the 2.87 meter results shown in Fig. 14. However, 
in both of these cases the bimodal character exhibited by the experimental 
data is well represented, with the spatial location (if not the magnitude) 
of the temperature peak agreeing reasonably well with that measured in the 
experiments. At a height of 5.73 meters, the peak temperature is still 
overpredicted, and the predicted bimodal structure is more pronounced than 
that shown by the experimental data, as can be seen from Fig. 15. At this 
vertical location the predicted temperature profile is also considerably 
wider than the data indicate. This anomaly was observed in all calculations 
of the 15 meter (50 foot) diameter fire that were carried out: although the 
temperature profile width is in general well represented for measurement 
stations other than that at 5.73 meters, it is not well represented at 5.73 
meters. One possible cause of this anomaly could be an offset in the experi- 
mental centerline at this point, but note that the data from both of the 
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15 meter (50 foot) fire tests overlay reasonably well at this point. Another 
possible cause may be the accumulated effects of the two rows of sample elements 
located below this height, but a mechanism that would produce the observed 
effects has not been identified. 

It should be recalled from the discussion in Section 2.8 that each of 
the profiles so far described is in the region in which the approximations 
inherent in the parabolic formulation can, on the basis of computed stream 
function patterns alone, be expected to be the least applicable. In fact, it 
was observed in Section 2.8, on the basis of preliminary computations in which 
higher predicted temperatures and thus larger effects of buoyancy than in the 
model currently being discussed were evident, that from the stream function 
pattern it did not appear that reasonable results would be achievable with a 
parabolic model at distances less than two or three pool radii from the 
liquid surface. However, the results shown in Figs. 16 and 17 indicate that 
very good radial profile results are, in fact, achieved for h * 11.48 and 21.43 

meters, or h/r = 1.54 and 2.86, respectively. 

P 

3.2 SPECIES PREDICTIONS AND COMPARISON WITH EXPERIMENT 

A variety of data, encompassing both radial and axial profiles, is available 
for the chemical species observed in the 15 meter (50 foot) diameter fire tests. 
The data include measurements of O 2 , N 2 , CO, CO 2 , soot and unburned fuel. 

It should, however, be noted at the outset of these comparisons that the 
chemistry model used in the pool fire code was, as described in Section 2.8, 
designed primarily to obtain reasonable temperature predictions, and not to 
predict in detail the species mass fraction distributions within the fire plume. 
Indeed, as will be seen, the experimental data indicate a considerably slower 
overall rate of reaction than that predicted by the quasi -equilibrium model 
utilized in the pool fire analysis. Nevertheless, a comparison of predicted 
and measured species distributions is of some interest. 

Detailed radial profile data are available for the 5.73 meter height, 
and the measured and predicted O 2 and N 2 profiles at this height are compared 
in Figs. 18 and 19. It will be recalled from the discussion of the temperature 
comparisons that at this height the thermal extent of the fire was overpredicted, 
and it can be seen from Figs. 18 and 19 that at this height the width 
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FIGURE 16. Radial Temperature Profiles at 11.5 M 
15 M (50 ft) Diameter Pool Fire 


of the fire plume as evidenced by both the O 2 and N 2 concentrations is over> 
predict d as well. Note that in the computations the turbulent Prandtl and 
Lewis numbers were both assumed to be unity, so that the element mass fraction 
and total enthalpy prefile widths will be equivalont arj identical to the 
velocity profile width as well. As was the case for the temperature profile, 
although the magnitude of these concentration profiles is not well predicted, 
the general shape of the profiles is accurate (di.sregarding the possibly 
anomalous centerline nitrogen mass fraction value shown on Fig. 19). 

Centerline values of the fuel mass fraction are of considerable interest 
in assessing the overall behavior of the model, and in interpreting the 
experimental results. A comparison of computed centerline fuel concentration 
wit!i experimental data for this parameter is shown on Fig. 20. The data 


48 


®***®*nal PAor n 
OF POOR QUAUiy 



TESTS 


TESTS 


^ FIGURE 17. Radial Temperature Profiles At 2.14 M 


IS M (50 ft) Diameter Pool Fire 


indicate persistence of fuel along the centerline until about two pool radii 
above the pool surface, and this long persistence of fuol is not predicted 
by the computations. Note again, from Fig. 20, the improvement in the predic- 
tion relative to the experimental data that is generated by the use of the 
unmixedness model. Although the unmixedness model used in these computations 
is by no means a rigorous treatment of the phenomenon, but is rather an 
order-of-magnitude estimate, the fact that even with the incorporation of 
unmixedness the computed fuel concentration is, below about 11 meters, consider- 
ably lower than that measured, provides evidence that the fuel consumption rate 
is much slower than would be predicted by equilibrium chemistry. 

It will be recalled that for temperatures above 750*K, a major branch 
of the quasi-complete-equilibrium chemistry model involves the production of 
CO and soot, at fuel -rich equivalence ratios. This branch of the model is 
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IE 18. Comparison of Predicted and Measured O 2 Radial 
Profiles, 15 M (50 ft) Diameter Fire 


FIRE DIAMETER 15m (50 ft) 



FIGURE 19. Comparison of Predicted and Measured N 2 Radial 
Profiles, 15 M (50 ft) Diameter Fire 
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/ FIGURE 20. Comparison of Predicted and Measured Centerline 
Fuel Concentration, 15 M (50 ft) Diameter Fire. 

reached relatively early in the fire development, as evidenced by the CO and 
soot mass fraction profiles shovm in Figs. 21 and 22. As the local equivalence 
ratio drops, the CO 2 -H 2 branch is activated, and as the equivalence ratio 
becomes less than unity, the products of combustion become CO 2 and H 2 O. The 
advent of these branches of the model is shown on Fig. 23. Both the early 
CO and soot mass fractions are strongly overpredicted, which is consistent 
with the rapid disappearance of fuel shown in Fig. 20. Taken together, these 
three comparisons provide a further indication that the fuel consumption in 
the fire is highly nonequilibrium process, as well as one which is strongly 
affected by turbulence-chemistry interaction or unmixedness phenomena. 

At the conclusion of this comparison, it should be reiterated that the 
chemistry model used in the fire model is intended to provide reasonable temper- 
ature estimates but not necessarily to predict the observed species distribu- 
tions. Indeed, the measured species data obtained in the IS meter (50 foot) 
fire experiments strongly indicate the nonequilibrium nature of the chemical 
processes occurring within the fire. A quasi-complete equilibrium formulation 
cannot be expected to adequately model the processes that are observed. 
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FIGURE 21. Comparison of Predicted and Measured Centerline 
CO Concentration, 15 M Fire. 

3.3 RESULTS FOR OTHER FIRE SIZES AND FIRE SCALING 

In addition to the computations of the 15 meter (50 ft) and fire, computa- 
tions were also carried out for a 7.6 meter (25 foot) and a 3; 8 meter (12.5 ft) 
diameter pool fire. In general, the results obtained for the smaller sizes 
were similar to those already discussed for the 15 meter (50 ft) diameter 
fire, and these results are summarized in the form of isotherm and iso-soot- 
concentration contour plots in Figs. 24-29. These figures show that the basic 
features of each of the fires are quite similar, with a characteristically 
bimodal temperature distribution being seen at the base of the fire, followed 
by a high- temperature plume. Qualitatively, the isotherm distribution that is 
computed using the parabolic model is also quite similar to that observed 
experimentally. Fig. 30, although there are some quantitative differences. 

In general, however, the computed results provide a good overall picture of 
the thermal environment within the fire. Soot distribution profiles are shown 
in Figs. 25, 27, and 29; although the comparisons with experimental data 
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FIGURE 22. Comparison of Predicted and Measured Soot Concentrations 
Along Fire Centerline, 15 M (50 ft) Diameter Fire 

described in the previous section indicate that the model overpredicts the 
initial soot concentration, the comparison of results for the different fire 
sizes indicates that the soot concentration profiles for each of the fire 
sizes are also qualitatively similar. Fig. 31 is an overlay of the computed 
isotherms for the 15 meter (50 ft) diameter fire on a drawing made from a 
photograph of one of the fire tests; since the detailed comparisons indicate 
that the computations are in general agreement with experimental results for 
this fire size, this overlay provides an interesting depiction of the flame 
structure. 

For many turbulent flows, geometric scaling holds provided the Reynolds 
number is large enough; i.e., for nozzles twice the diameter, the length to 
a given level of mixing is twice as long. Rigorously, this scaling applies 
only when the source terms in the governing equations are either zero or 
negligibly small. This is not the case for these fire calculations, in which 
the buoyancy terra provides a nonnegligible source term for the momentum 
equations. However, as Fig. 32 shows, the centerline temperature profiles for 
the 15 meter and 7.6 meter fires scale reasonably well, while results for the 
smallest fire size studied fall below those for the other two sizes. Note that 
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FIGURE 23. Comparison of Predicted and Measured Centerline 
CO 2 Concentrations, IS M Fire. 

the axial location of the peak centerline temperature scales almost exactly, 

although it should be recalled that in the region of peak centerline temperature 

% 

the assumptions basic to the use of a parabolic fire plume model are not well 
met . 


Centerline velocity profiles are also of interest, both from the standpoint 
of scaling and because the magnitudes of the velocities encountered in fire 
plumes of this large scale often tend to be underestimated. Fig, 33 shows 
centerline velocity results for the three fire sizes. Here it is evident that 
geometric scaling does not apply, and that velocities of the order of 30 m/sec 
can be encountered in large-scale fires. The data point indicates on Fig. 33 
was obtained during the NASA/White Sands test program, for a 15 m (50 ft) 
diameter fire, using an anemometer designed by A. D. Little, Inc. to function 
under extreme ambient conditions existing within the fire. The predicted 
centerline velocity for the 15 m (50 ft) fire is somewhat higher than this 
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FIGURE 24. Predicted Isotherm Contours, "K, For 3.8 M 
(12.5 ft) Diameter Pool Fire 
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/ FIGURE 28, Predicted Isotherm Contours, (*K), 15.2 M 
(50 ft) Diameter Pool Fire 


measurement indicates, but in the absence of an error estimate for the 
instrument the level of agreement shown can be considered to be satisfactory 
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FIGURE 29. Predicted Soot Mass Fraction Contours 
15. 2 M C50 ft) Diameter Pool Fire 


60 


HEIGHT 





soft FIRE TESTS 


24 r 



2 4 6 

RADIUS, m 


TEMP.,'K 

cggg >1500 

1400-1500 
fTTim 1300-1400 
C=za 1200-1300 
rrm. 1100-1200 
1000-1100 
800-1000 
I — ^ 400-800 


j 

8 


FIGURE 30. Isotherm Contours Obtained From Measurements 
in 15.24 M (50 ft) Fire 
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FIGURE 32. Centerline Temperature as a Function 
of Scaled Distance From Pool 
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FIGURE 33. Centerline Velocity as a Function of 
Scaled Distance From Pool 
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4. PREDICTION OF FIBER FLUXES AND CONSUMPTION 


The model described in the preceding section has been applied to the 
computation of carbon fiber transport and consumption within a fire plume. 
Similar initial conditions were used for all fire calculations: fuel vaporiza- 

_3 

tion rate of 8.5 x 10 cm/sec (0.20 in/min), particle release rate of 
10 kg/min, distributed uniformly across the pool, a fuel vapor temperature of 
400®K and an ambient temperature of 300®K. The particle characteristics 
for these computations were as follows: 


Particle Class 

I 

rp(cm) 

tp(cm) 

fp (kg/min) 

Remarks ^ 

1 

3.5x10"^ 

0.15 

0.288 

Single Fiber 

2 

3.5x10"'* 

0.35 

0.448 

Single Fiber 

3 

3.5x10"* 

0.65 

0.791 

Single Fiber 

4 

3.5x10"* 

1.00 

1.217 

Single Fiber 

5 

3.5x10"* 

1.40 

1.076 

Single Fiber 

6 

3.5x10'* 

5.00 

3.202 

Single Fiber 

7 

3.5x10"* 

25.00 

1.282 

Open Fiber Cluster 
Long Single Fiber 

8 

5.12x10"^ 

10.00 

1.694 

Large Fiber Cluster 


4.1 FIBER CONSUMPTION PREDICTIONS 

Fiber consumption rates were obtained as part of the work described in 
Ref. 4; although this parameter was not computed in more recent work, a 
comparison of the fire structure results presented herein with those presented 
in Ref. 4 indicates that the modifications made to the chemistry model in this 
work would not materially affect the fiber consumption results. In view 
of the thermal nonequilibrium results reported in Section 2.7, the thermal 
equilibriiim assumption was retained for these calculations. Results, expressed 
as percent of total fibers consumed, are shown for the single fibers in 
Fig. 34: single fiber results are shown because these fibers are much more 

readily characterized in terms of geometric parameters (length, diameter, 
effective surface area, volume) than are the fiber clusters and clumps. These 
results show that the bulk of the fiber consumption occurs in the lower regions 
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FIGURE 34. Fiber Consumption (Single Fibers) Computed 
For Two Fire Sizes 

of the fire. For the larger fire, the initial consumption at a given value 
of height above the liquid fuel pool is lower than for the smaller fire since 
the height at which the peak temperature occurs increases with pool diameter. 
It is also evident that the total fiber consumption increases somewhat with 
fire size, and that in both cases fiber consumption has ceased by the 60 meter 
(200 ft) height in the fire. 

4.2 FIBER FLUX PREDICTIONS 

Fiber flux results for the single fibers are presented /or three fire 
sizes in Figs. 35-37. These results are presented as local fiber fluxes 
(kg/m -sec) normalized with initial fiber flux: since the burning rate for 

all single fibers is the same, because of the assumption that the surface 
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area of the ends of a cylindrical fiber can be neglected in the consumption 
calculation, and it has been assumed in these computations that all single 
fibers have the same turbulent diffusivity, these normalized curves can be 
used to relate the flux of any class of single fibers to its initial flux. 

The initial fiber flow rates used to generate these computational results 
have been tabulated at the beginning of this section. 

In all cases, because of the tendency of the fire to "neck down," the 
large majority of fibers are observed, even at the 61 meter (200 foot) height, 
to be in a region within one pool radius of the centerline. In addition, 
this "necking down" tendency, as shown in Section 2.5.1, leads to high-altitude 
( — 100 ft or greater) fiber concentration profiles that are essentially 
independent of the initial location of the fibers within the fuel pool at 
h = 0; therefore, the results shown in Figs. 35-37 can be used to estimate 
fiber fluxes for any arbitrary (but axisymraetric) initial fiber distribution. 
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FIGURE 35. Normalized Fiber Flux Profiles as a Function of Height, 
3.81 M (12.5 foot) Diameter Fire 
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FIGURE 37. Normalized Fiber Flux Profiles as a Function 
of Height, 15 M (50 foot) Diameter Fire 


5. CONCLUSIONS 


The result of the work outlined in this report has been the development 
of a model of pool fire phenomena which, although approximate, is capable of 
providing considerable detail of the structure of a pool fire in zero-wind 
conditions, particularly with respect to expected temperature distributions. 
Because of the approximations inherent in a parabolic formulation of the 
pool fire problem, there are serious discrepancies between the model predictions 
and experimental data at low altitudes in the fire plume, but the level of 
accuracy of the model results increases markedly with vertical distance from 
the ground plane. While there is only a limited amount of experimental data 

i 

available with which to compare model predictions, these data provide the 
most complete experimental data on large-scale pool fires that is currently 
available. 

Specific results of the work described in this report are as follows: 

1. An analytical model of large pool fires which is in good agreement 
with measured fire structure data has been developed. This model 
can be used as a tool for the analysis and interpretation of large- 
scale fire measurements and consumption of fibers in a large-scale 
fire. 

2. Measured and computed temperature and species profiles show that 
significant unmixedness and incomplete combustion effects are present 
in large JP-4 pool fires. 

3. Computational results indicate that the effects on fiber consumption 
of locally high temperatures because of temperature fluctuations 

in the fire are negligible; similarly the short residence time of 
fibers in regions of the fire calculation in which mean t«nperature 
is overpredicted compared to the measurements indicates that the 
overprediction of temperature in these regions should not greatly 
affect the overall fiber consumption prediction. There is, however, 
a need for further verification of the fiber consumption rate expres- 
sion over a range of oxygen concentrations typical of a fire. 


4. While particle thermal nonequilibrium can result in some tempera- 
ture mismatch between the fibers and gas phase* the overall effects 
of nonequilibrium particle heat-up and cool-down tend to compensate 
resulting in only a small overall effect on particle consumption. 

5. Based on the particle consumption expression defined in this report* 
overall single-fiber particle consumption of about 15% of the total 
particle mass is computed for large pool fires. 

Further improvement of the mathematical modeling of the structure of a 
large pool fire over that described in this document requires recourse to a 
two- or three-dimensional elliptic Navier-Stokes formulation: beyond the 

complexity that such a formulation introduces of itself* establishment of 
appropriate boundary conditions for a more exact model is a difficult task. 
However. The use of an elliptic approach will allow the relaxation of the 
assumptions related to mixing rate that are incorporated in the present 
formulation* and such a model is the only valid method of treating in detail 
the phenomena in the lower part of the fire. 
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